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Abstract 

In this note radiative corrections to the total hadronic decay rate of 
the r-lepton are studied employing perturbative QCD and the opera- 
tor product expansion. We calculate quadratic quark mass corrections 
to the decay rate ration to the order 0{a1rn?) and find that they 
contribute appreciably to the Cabbibo supressed decay modes of the 
r-lepton. We also discuss corrections of mass dimension D = 4, where 
we emphasize the need of a suitable choice of the renormalization scale 
of the quark and gluon condensates. 
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1 Introduction 



Within the three lepton generations which to our present knowledge constitute the 
leptonic sector of the Standard Model, the r-lepton is the only particle decaying 
into a semihadronic final state. As was pointed out ^, ^, ^ |, |[ some 
time ago, the methods of perturbative QCD can be applied to estimate the decay 
rate ratio 

r(r ^ i/^hadrons) 

r(r UrCVe) 

Work has also concentrated on nonperturbative effects 0, |, |^, |[ as well as on 
electroweak corrections ^ to this quantity. 

In order to calculate QCD corrections to Rr we consider the 2-point correlators 
for the vector (j^j^ = qi'y^qj) and the axial vector (j^^^- = qi'y^'y^qj) currents 
= u,d,s): 



n;^,,(g) = ^ J dxe^'^^(j\j;!-ix)j:^pm) (2) 

As usual the correlation functions may be decomposed into a transversal and a 
longitudinal part 

nj'/^(g) = {-g,uq' + q,qu)n^yliQ') + q,qunP/^iq'), (3) 

where the spectral functions nyj^(g^), nyj^(g^) correspond to hadronic final states 
with respective angular momenta J = 0, J = 1 in the hadronic rest frame. The 
hadronic decay rate of the r-lepton is obtained by integrating the absorptive parts 
of the spectral functions with respect to the invariant hadronic mass: 



ds S 



1 + 2— j Imn(^) (s + ie) + ImH^^^ (s + ie) 



(4) 
where 

n(^) = \Vu,nii% + nlZ) + \Vusnii[i]v + ng^). (5) 

Due to the analyticity of the spectral functions in the complex s-plane, which is 
cut along the real positive s-axis, Rr can be expressed as the contour integral 
along a circle C of radius |s| = M^: 

2 



J\s\=M? M2 I M2 i 



1 + 2^ n(-^)(.)-2^n(o)(.) 



(6) 



l\s\=M?M^\ M^J 

We have used in this equation the combination 11'^''"''^^ (g^) = Il^^\q'^) + ]l^^\q'^). 

Rr may be expressed as the sum of different contributions corresponding to 
Cabbibo suppressed or allowed decay modes, vector or axial vector contributions 
and the mass dimension of the corrections: 

Rr = Rr,v + Rt,a + Rt,s (7) 



1 



with 
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E 






D=2,4, 
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^ D=2,4,. 



'A,ud I ) 



(8) 



Rs= 3|K.r 1 + 5^°^+ E 



Here D indicates the mass dimension of the fractional corrections ^\fv/A ^^"^ 
dlf"* denotes the average of the vector and the axial vector contributions: dlf' = 

In the literature results for perturbative QCD corrections to current-current 
correlators have been published by several groups. For the electromagnetic current 
correlator corrections were calculated in the massless limit up to 0{a'^) [|ri|] and 



0{al) [jT2|, whereas quadratic mass corrections to this quantity were given to 
second Ola"^) [0 and third 0{al) order fl^. For the axial current correlator 



mass corrections are known from two loop |]T5|, [16|, [T^] and three loop calculations 



IS] for flavour non singlet type diagrams. Flavour singlet contributions to the 



axial correlator involving purely gluonic intermediate states were studied in the 



heavy top limit in |T9|, Q and for the massive case in |21[]. Finally second order 
massless corrections for the scalar current correlator can be found in ||22|| . 

A detailed analysis of QCD corrections for the r decay into hadrons has been 
performed in where much attention focussed on nonperturbative contributions 
to corrections of higher dimension. 

The aim of this letter is twofold. First we extend the analysis of ref. ^ and 
calculate the order C(a^) term of the (dimension D = 2) longitudinal spectral 
function nyj^(g^) in its power expansion with respect to as- This analysis is done 
in section 2. 

The second part of this note is contained in section 3 where we consider corrections 
of mass dimension D = 4 to Rr- We estimate their size where we use a scale fi 
for the quark condensates which corresponds to the energy scale of the process 
under consideration. Our choice /i^ = seems to be more appropriate to us 
than /i^ = oo as used in p. 

Our numerical results are discussed in section 4. We finally list some formulae in 
the appendix. 
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2 Dimension D = 2 Corrections 



The way to derive the quadratic mass correction to the longitudinal spectral 
function Tlyj^.-j is based on the knowledge of the anomalous dimensions 7^^'''^^ 



of the vector and the axial vector correlators as defined in ||T8 

,2 



(9) 



Up to order 0{a1) the anomalous dimension 7^^'''^^ is given by (with / denoting 
the number of quark flavours) 



7 



r/^-^ = 6 1 + 



5a^(/i^) /as(/i^) 



3 TX 



+ 



TT 



(10) 



and governs the renormalization group equation for the longitudinal spectral func- 
tion 



d 



r{0) 



^ {rriiifi ) ^mj(/i )) w/aa 1 

2 ■^^y/A;jj Q2 !■ 



(0) 



IGtt^ 



In order to get ^yj^ ^^e RGE eq.(||)must be integrated. The solution reads 



n 



(0) 

V/A;ij 



1 {mi{fi^)Tmj{fi^)Y 



exp 



dA' ^" 



VV/AA 



exp 



_2 / , dA"^ 



:i2) 



where the /^-function and the quark mass anomalous dimension 7^ are given in 
the appendix. The integration constant C(/io) may be fixed by a specific choice 
for the arbitrary scale /Xg which we choose to be fi^ = Q"^. This leads to 



TT 



2\ r 1Q1 

4C(3) - f 



(13) 



With the relation between the //^-dependent mass of the MS-scheme and the 
Q^-dependent running mass (see appendix) 



m (/i ) = m (/Iq) exp 



-2 dA 



7m 



(14) 



we cast Uyj^.-j into the following form 



n 



(0) ^ 1 jm^mTm.my 

V/A;ij ^g^2 Q2 



exp 



dA 



/1m 

T 



dA'- 



VV/AA 



1 {mi{fi^)Tmj{^^)Y 
167r2 Q2 



exp 



exp 



-2 



A' 



dA"^ 



-C{Q' 



dA 



(15) 



as(l£ 



ryVV/AA 

dA exp 

P 



dA'^ 



Inserting the series expansion for /? and 7^, as given in the appendix we obtain 
after integration 



n 



(0) ^ 

V/A;ij 271-2 



3 K(Q2)Tm,(Q2))2 a,{Q 



5 
2 



+ 



''''' + SC(3) 



4896 34 
3 (mi(/i2) ipmj(/i2))2 J /Q;^(yri2^\ -1 



TC 



TT 



+ 



8011 41 
'4896 ^34^*- ! 



TT 



(16) 



The leading and next to leading terms of this expansion agree with p whereas 
the coefficient of as is new. The /i^-dependent part of eq. (|16|) does not contribute 
to the hadronic r-decay rate. 

Mass corrections of the transversal spectral function n^^^^]^. were calculated 
to order 0{a1m?) for the electromagnetic correlation function in |jl3|. Recently 
one of us (K. Ch.) extended this calculation to the case of unequal masses. The 



A 



obtained result reads 



^W/A;ij 



3 {m,{Q^)±m,{Q^)f 



87r2 



1 + 



+ 



1ZH51 + 5Hc(3) - ^c(5) 

432 27^^ ' 54 , 



TT 



3 K(Q2)^^^.(g2))^ 



87r2 



1 + 2 



+ 



vr 

4351 13 



144 



1 1 5 

3C(3)--C(5) 



TC 



127r2 



TT 



(32 - 24C(3))E 



(17) 

Even for the decay of the r-lepton into non strange quarks is a dependence on 
rris introduced by the second order contributions due to the s-quark circulating 
in a fermion loop. The term proportional to J2k ^kiQ"^) / was first computed in 

nij the result eq. ([T7|) exactly reproduces 



24] For the vector correlator with rrii 
the one of [|13|. This means that a slightly different result (namely for the terms 
proportional to C(3)) for the diagonal vector correlator given in [^] seems to be 



wrong. For a detailed discussion we refer the reader to |23 



3 Dimension D = 4 Corrections 



Besides the perturbative radiative corrections to Rr also nonperturbative QCD 
effects influence the hadronic r decay rate. The short distance operator product 
expansion (OPE) for the spectral functions 

1 



n(^)(- 



E 

D=0,2,4,. 



)D/2 



1^ 



dimO=D 



may be used to take into account both perturbative and non perturbative contri- 
butions. For the fractional corrections 51^1^ (^^^ 6q.(P)) we obtain 



E 



£111110 = ^ 



MP 



ds f —s 



-D/2 



S 



(19) 

The local operators in this expansion for the D = 4 perturbative corrections are 
the unit operator multiplied by a quartic product of quark masses and the vac- 
uum expectation values of composite operators constructed with gluon and quark 



fields: {GG{fi)), {mi^jj)^ . The latter contain nonperturbative contribu- 
tions p6|, as well as mass logarithms of the form ||T6| , ^ m^Q;" (/i) ln^(m//i) 
and depend nontrivially on the renormalization scale jU via the corresponding 



renormalization group equations [g^. As usual, setting the renormalization scale 
/i = Q in eq.(p!8|) allows to absorb all logarithms In /i^/Q^ appearing in the coeffi- 
cient functions C'^'^-* into the running coupling constant. This procedure also leads 
to an implicit Q^-dependence of VEV's, which is not convenient for a numerical 
analysis. The common remedy is to to solve the corresponding RG equations and 
express {O) in terms of as{Q) and some RG invariant combination not depending 
on Q2. 



As has been shown in [28| it is possible to construct linear combinations of 



the operators which are invariant with respect to an arbitrary scale fi: 
da) ^(l + }^^]^(GG(^)) 



371 



These combinations were intensively used in the analysis of |^ under the names of 
the gluon condensate {^GG) and the quark condensate (mj^j\l/j) respectively. 
Of course the choice is by no means unique. Moreover, it directly leads to an 
l/os enhancement factor in the dimension D = 4 correction to the ratio Rr 
(see eq.(3.10) in P). This results in a partial cancellation between quartic mass 
corrections and those coming from the quark condensates, which in the final 
analysis causes some partial loss of accuracy of theoretical predictions 0). 
In addition one can easily check that 

{Ig)= lim^(GG(A)), 
{Xi,)= hm(m,(^)'I',-^,(/i)). 



Thus the choice of eq.(POD as RG invariant vacuum condensates features a quite 
high normalization scale for operators involved in a problem with a typical mo- 
mentum transfer of about 1 GeV. No wonder that this causes an artificially large 



perturbative mass correction of order 0{m'^) to various 2-point correlators as was 
found in |15| . 

Keeping all this in mind we have chosen just the very quark and gluon conden- 
sates normalized at the "natural" scale fi = Mr as our RG invariant condensates: 



TT 



GG) 



TT 



(22) 



(171- 



Using the RG invariance property of the combinations eq.(pO|) it is a straightfor- 
ward exercise to express the Q^-dependent VEV's in terms of our RG invariants 
and the running as{Q),m{Q). 

The result for the fractional corrections after performing the contour integral 
reads 



r(D=4) _ ^4 



11 



-TT 



TT 



TT 



GG) 



-IStt' 



9 fas{Mr 



TT 



TT 



84 



'as{Mr) 



TT 



(23) 



- — [mi{Mr)Tmj{Mr)][mi{Mr)Tm%Mr)] 
±Qmi{Mr)mj{Mr)[mi{Mr) =F mj{Mr)f 



+mml{Mr)m]{Mr) 

For some terms the leading coefficient is of order 0{a1) due to the fact that a 
term of given order of the transversal spectral function is contributing only in 
higher orders to the fractional correction and order (9(a°) terms may eventually 
integrate to zero. We have neglected terms of order a^m^ in eq.(P3D. No enhanced 
term proportional to the inverse power of the coupling constant occurs as was the 
case for the choice of reference 0. 



4 Numerical Discussion 

The numerical discussion in this section is based on the parameters (quark masses, 
condensates etc.) as they are given in the appendix. Depending on the value 



7 



of Aqcd the running coupling constant as{M^) ranges between 0.16 and 0.44 
and is thus sufficiently small for the perturbation expansion to be meaning- 
ful. The quadratic mass corrections (^[^Xi] obtained from the contour integra- 
tion of the spectral functions eqs.(|l^,0) are collected in Table One observes 
that corrections from nonstrange decays are negligible, whereas the mass of the 
strange quark affects the ratio Rs for strange decays by —20% for an intermediate 
as{M^) = 0.3. When the total hadronic ratio Rr is considered, this rather large 
strange mass contribution is reduced by the Cabbibo supression factor iKisP to 
—0.8% {Vud = 0.9753, Vus = 0.221). The influence of the second order correc- 
tion depends of course on the value of Aqcd- Even for nonstrange decays strange 
quark mass effects are present at order 0{a'^) due to a virtual strange quark 
loop. They are quite comparable in size to the leading mass corrections of order 
(m^im^)^ and enter with opposite sign. Due to their increasing size for larger 
the corrections ^Iya w show little dependence on Aqcd- For strange decays of the 
r-lepton we recall the numerical values of the coefficients entering the spectral 
functions (in the limit of vanishing masses of the light quarks) 



(0) _ 3 m|(Q^//a^V'_5_ a,(Q 

^W/A;us- 2772 Q2 I TT ] 2 ^'^^^ TT 




2^ 



3 m'',{fi^) ] f as{fi^)\ ' 2 - 1 867-"'^^'^ 



(24) 



27r2 Q2 I / ■ TT 



2' 



and the fractional correction 



2 



- -8^^ \l^'4"-^^ 46.002 f n . (26) 



TT \ TT 




Second order contributions to 5(^=2) originate not only from order 0(«s) terms 
of the spectral functions, but are also induced by lower order terms. 

In Table 2 the corrections of mass dimension D = 4 are shown. The change 
of scale from fi = 00 to fi = Mj. significantly affects only the Cabbibo supressed 
vector contribution. The small corrections Sy^us^^ in the analysis of were due 
to a numerical cancellation between the leading quark condensate and an l/a^- 
enhanced mass term m^. This resulted in a large relative uncertainty for this 
correction. When the condensates are defined at the scale fi = Mj. a similar 
compensation of numerically large terms does not occur due to the absence of 
the l/a^-enhanced mass term. Compared to the corresponding numbers of ^ 



We thank K. Maltman |29 for pointing out an error in the original version of this paper. 



the fractional corrections 6y'.~g^ are increased by an order of magnitude which 
reduces the relative uncertainty considerably. In addition they do not change their 
sign at large values of Aqcd- The strange D = 4 corrections contribute only little 
(namely -3%) to the strange decay ratio Rs and are negligible when the total 
r-decay ratio Rr is considered. 

We now include all corrections to present the values for the separate contribu- 
tions as well as for the total hadronic r-decay rate in Table 3. Here electroweak as 
well as corrections of higher mass dimension have been taken into account from 
P, where it has been pointed out that the biggest nonperturbative corrections 
to Rr in fact come from dimension D = 6 condensates. Our estimation of the 
corresponding uncertainty on R^- is more conservative than in because the 
large cancellation between the vector and the axial vector contribution is based 
on the validity of vacuum saturation approximation and may not necessarily be 
translated into a similar compensation for the uncertainties of the separate con- 
tributions. Numerically Rr is dominated by purely perturbative contributions of 
zero mass dimensions which survive in the massless limit. According to ||3^ a 
proper summation of the effects of analytical continuation from space-like mo- 
menta to time-like ones is a necessity for these terms. At the moment it is not 
quite clear for us how essential these effects are for power suppressed contribu- 
tions that we are dealing with. We hope to return to this and related problems 
in future publications. 

To conclude, we have studied power suppressed perturbative and nonpertur- 
bative QCD corrections to the semileptonic decay rate of the r-lepton. Strange 
quark mass effects contribute considerably to the decay modes with strangeness 
content and still reach the percent level for the total hadronic decay rate. Sec- 
ond order 0{a'^) corrections introduce a strange mass dependence even for non 
strange decays. We also have studied corrections of mass dimension D = 4, where 
we discussed an appropriate choice for the renormalization scale of the quark and 
gluon condensates. In view of the accuracy of QCD predictions for the decay rate 
of the r-lepton into hadrons, semileptonic r decays remain an important and 
interesting tool for testing QCD. 



Q 



Table 1: S fractional corrections 





as{Mr) 


dD=2) _ . .3 
°ud,V -"-^ 


• 10^ 


dD=2) 


r(0=2) 


52 


0.16 


-0.20 ± 0.04 


-0.37 ± 0.07 


-0.074 ± 0.017 


-0.077± 0.017 


308 


0.30 


-0.36 ± 0.12 


-0.87 ± 0.20 


-0.197 ± 0.044 


-0.206 ± 0.045 


547 


0.44 


-0.41 ± 0.29 


-1.60 ± 0.46 


-0.418 ± 0.094 


-0.439 ± 0.096 



Table 2: 5--TriA fractional corrections 



AgcD/MeV 












52 


0.16 


0.17 ± 0.13 


-5.00 ± 0.61 


0.011 ± 0.005 


-0.048 ± 0.007 


308 


0.30 


0.60 ± 0.27 


-4.68 ± 0.67 


0.010 ± 0.005 


-0.050 ± 0.008 


547 


0.44 


1.30 ± 0.55 


- 4.20 ± 0.83 


0.009 ± 0.005 


-0.055 ± 0.009 



Table 3: Contributions to the hadronic r decay rate Rr 



AgcD/MeV 


as{Mr) 


Rv 


Ra 


Rs 


Rt 


52 


0.16 


1.59 ± 0.02 


1.49 ± 0.03 


0.145 ± 0.003 


3.23 ± 0.02 


308 


0.30 


1.73 ± 0.02 


1.63 ± 0.03 


0.140 ± 0.007 


3.51 ± 0.05 


547 


0.44 


1.95 ± 0.08 


1.85 ± 0.08 


0.128 ± 0.016 


3.93 ± 0.16 
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A Appendix 



The /9-function and the quark mass anomalous dimension 7^ are defined in the 
usual way 

(27) 



<l^^ i>o 



a. 



TX 



(28) 



Their expansion coefficients up to order 0{aj.) are well known |32| and read 



/3o =(ll-^/)/4, /3i = (^102 - /16, 



(29) 



7o = 1, 7i 



202 20 



72 



( 1249 



2216 160,, , 



(30) 



The Riemann zeta function has the values C(3) = 1.2020569,^(5) = 1.036927. 
The solution of eq. (^Tl) is given by (L = ln/i^/A|^) 



as(/i^ 



TT 



■ 1 AlnL 1 



(31) 



while eq. (EBp is solved by 



m(/i^) = rh 2(3q 



2\\ 70 /A) 



1 + 



as(/i^ 



TT 



TT 



7l _ Pilo 



+ 



1 / a,(/i2)y 



TT 



7i _ A70 
/5o /5o' 



_^72 _ /?i7i _ /?27o _^ l^ho 
Po Po Po Po 



(32) 



For the renormalization invariant quark mass parameters we have taken the values 



33, 34 



m„ = (8.7± 1.5)MeV", = (15.4 ± 1.5)MeV, m, = (270. ± 30)MeV. (33) 
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As input value for the gluon condensate we have used 



{asGG) = (0.02 ± OM)GeV\ 



(34) 



whereas the quark mass condensates are parametrized by |Q (mj\I/j\&j) = —rhifi^ 
with 
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